Fragility in p-spin models 
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We investigate the relation between fragility and phase space properties - such as the distribution 
of states - in the mean field p-spin model, a solvable model that has been frequently used in studies of 
the glass transition. By direct computation of all the relevant quantities, we find that: i) the recently 
observed correlation between fragility and vibrational properties at low temperature is present in 
this model; li) the total number of states is a decreasing function of fragility, at variance of what 
is currently believed. We explain these findings by taking into account the contribution to fragility 
coming from the transition paths between different states. Finally, we propose a geometric picture of 
the phase space that explains the correlation between properties of the transition paths, distribution 
of states and their vibrational properties. However, our analysis may not apply to strong systems 
where inflection points in the configurational entropy as a function of the temperature are found. 

PACS numbers: 64.70.Pf, 61.43. Fs, 75.50.Lk 



I. INTRODUCTION 



The glass-forming materials are characterized by a 
huge variation of their transport properties (viscosity, 
mobility, diffusivity) and dynamical (relaxation times) 
properties, upon supercooling. As an example, the vis- 
cosity in the whole "liquid" range, if crystallization is 
avoided, increases of about 17 orders of magnitude by de- 
creasing temperature before the system falls in the glassy 
state. The transition to the latter is conventionally fixed 
at the temperature, T g , where the viscosity reaches a 
value of 10 13 poise. Different systems show different tem- 
perature behavior of the viscosity, and they have been 
classified accordingly. Fragility is an index measuring the 
steepness of the viscosity as a function of the tempera- 
ture on approaching the glassy state: "fragile" systems 
(high value of the fragility index) are characterized by a 
super-Arrhenius behavior of the viscosity, that increase 
very fast and -if extrapolated below T g - seems to diverge 
at a finite temperature Tk- In "strong" systems (low 
value of the fragility index) , on the contrary, the viscos- 
ity increase is less dramatic and follows an Arrhenius law, 
apparently diverging only at zero temperature. 

The identification of the microscopic details that, in 
a given glass former, determine the temperature depen- 
dence of the viscosity, and thus the value of the fragility, 
is a long standing issue in the physics of supercooled liq- 
uids and glassy state. Large numerical and theoretical 
effort have been devoted to the attempt to relate the 
fragility to the specific interparticle interactions (e. g. 
strong glasses are often characterized by highly direc- 
tional covalent bonds, while the fragile one have more 
or less isotropic interactions). More recently, the atten- 
tion has been focused on the possible relation existing 



between the fragility and the features of the Potential 
Energy Landscape (PEL), more specifically the energy 
distribution of the minima of the PEL and the proper- 
ties of the basin of attractions of such a minima. With 
this respect, a key point is the existence of a relation 
between viscosity (or relaxation times) and the configu- 
rational entropy S(T) (i. e., the number of basins pop- 
ulated at given temperature), namely the Adam-Gibbs 
relation: 



V(T) = Voo exp 



\TE(T) 



(1) 



This relation has been extensively tested against the ex- 
perimental result, and it is now commonly accepted as 
" correct" . Despite its success in describing both numeri- 
cal and experimental data, the Adam-Gibbs relation has 
not been still derived in a clear way from microscopic 
models. This leaves unsolved the question of the micro- 
scopic interpretation of the parameter £ , that is usually 
believed to be related to the properties of the transi- 
tion paths between different minima of the potential en- 
ergy, such as the height of the barriers or the connectivity 
of the minima. By using the Adam-Gibbs relation, one 
could expect to relate fragility to the properties of S(T), 
i.e., to the distribution of basins in the phase space of 
the system. However, this possibility is frustrated by the 
lack of knowledge on the parameter £. Indeed, once a 
model for £(T) has been chosen, one can obtain the whole 
range of experimentally observed fragilities by varying £ 
(see Ref. [jj and references therein). More specifically, 
in Ref. Q it was observed that for a large class of mod- 
els for E(T) - where S(T) is a concave function of T 
that vanishes at a given temperature Tk and assumes its 
maximum E* at high temperature ( "Gaussian- like mod- 
els") - the relevant parameter that actually determines 
the fragility is 
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Thus, fragility appears to be determined by the ratio 
between £ (measured in units of fc^TV) and the total 
number of states E*/&b; it is related to both the distri- 
bution of minima (through £*) and the characteristic of 
the transition path between them (through £). The rela- 
tion between fragility and phase space properties can be 
even more complicated, in those cases where the function 
S(T) does not belong to the Gaussian class. 

The relevance of the concept of fragility also relies on 
the correlations that have been found between this in- 
dex and other properties of glass-forming liquids. Exam- 
ples of these correlations are the specific heat jump at 
T g (thermodynamic fragility) 0], the degree of stretch- 
ing in the non-exponential decay of the correlation func- 
tions in the liquid close to T g the visibility of the 
Boson peak at the glass transition temperature Q, or 
the temperature behavior of the shear elastic modulus in 
the supercooled liquid state [I| . More recently a striking 
correlation between fragility and the vibrational proper- 
ties of the glass at low temperatures has been found . 
Specifically, by examining the dynamic structure factors 
of different glass former well below T g , it has been found 
that the fragility of the corresponding liquid is propor- 
tional to the rate of change of the non-ergodicity factor 
in the T — > limit. Being the latter quantity fully deter- 
mined by the (harmonic) vibrational properties (eigen- 
modes of the disordered structure), this finding implies 
the existence of a deep relation between three features of 
the PEL: the energy of the minima, the transition paths 
between them (that together determine the fragility) and 
the Hessian matrix, evaluated at the minima themselves, 
that fixes the vibrational properties. 

With the aim to elucidate the existence of this unex- 
pected correlation between energy, curvature and tran- 
sition paths in the minima of the PEL we selected a 
solvable model of "glass", where i) the distribution of 
minima is "Gaussian-like" , ii) the vibrational properties 
of the minima can be determined, and iii) the transi- 
tion path between different minima can be evaluated and 
characterized by an energy parameter. More specifically, 
we investigate the mean field p-spin model (in both its 
spherical and fsing spin version) , a model that share with 
the structural glasses many aspects of the glass transition 
phenomenology, and that is known to have a Gaussian- 
like distribution of states. Our goal is twofold: i) we aim 
to verify if the analysis reported in Q is indeed correct 
in some microscopic model, i.e., if one can obtain a wide 
range of fragilities in a Gaussian-like model by varying 
the parameter £, and ii) to check whether one can explain 
the correlation between fragility of the liquid and the vi- 
brational properties of its glass found in 6] by studying 
the geometry of the phase space. The latter point could 
allow us to shed light on the origin of the correlation 
between number of minima, their vibrational properties 
and the property of the transition path between them. 

The paper is organized as follows: in Sec. [H] we de- 
fine the relevant quantities in the case of the mean field 
p-spin model; in Sec. lIIII we compute them for the spher- 



ical p-spin model, and in Sec. IIVI for the Ising p-spin 
model; in Sec.|3we discuss the relation between fragility 
and phase space geometry in these models, and compare 
our result with experimental data. Finally, we draw the 
conclusions. 



II. DEFINITION OF THE RELEVANT 
OBSERVABLES 

The quantities we wish to compute are (we will set kg — 
1 in the following): 

Tk Thermodynamical transition temperature 

T g Glass transition temperature 

Td Dynamical transition temperature 

E(T g ) Complexity at T g 

m(T g ) Fragility 

a(T g ) "Volume" of the equilibrium states at T g 

£{T g ) "Barrier height" at T g 

First, we have to identify the proper definition of these 
quantities in a mean field model. The main problem is 
that in a mean field model the glass transition tempera- 
ture T g is not a well defined quantity. Indeed, the relax- 
ation time of the system is known to diverge - as a power- 
law - when the temperature approaches the dynami- 
cal transition temperature Td, that corresponds to the 
usual Mode-Coupling temperature Tmct- The crossover 
from a power-law behavior of the relaxation time to an 
Arrhenius-like behavior, observed around Tmct in finite 
dimensional systems, is due to the activated processes 
becoming relevant; these processes are absent in mean 
field systems, and the crossover at Tmct becomes a true 
dynamical transition at Td \J\. To overcome this prob- 
lem, we will give an estimate of the height £{T) of the 
barrier that the system must pass through in order to 
escape from a metastable state at a given temperature 
T. Thus, we will make use of a "fictitious" Adam-Gibbs 
relation, 



V(T) = rjoc exp 



gen 

TE(T) 



(3) 



and define T g by r](Tg)/i] oc = const, or, equivalently, by 



T 9 Z(Tg 



= const 



(4) 



Note that in this paper we will not distinguish between 
the "complexity" (or "configurational entropy") E(T), 
that can be calculated in mean field models, and the "ex- 
cess entropy" measured in the experiments: indeed, the 
behave in a similar way in a wide class of systems 
Obviously, the quantity rj(T) has no dynamical meaning 
in a mean field context, but it provides an useful defini- 
tion of T g that hopefully coincides with the usual one in 
finite dimension. It will turn out that our analysis is not 
strictly dependent on this definition of T g , the behavior 
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FIG. 1: The two-replica potential for T€[Tg, T<j] in the spher- 
ical p-spin model. 



of the various quantities at T g being representative, as we 
will see, of a general trend observed at all temperatures 
T e [T K , T d ) by varying p. 



A. Two-replica potential 

The easiest way to provide a definition of all the above 
quantities is to consider the two-re plic a potential ap- 
proach introduced and discussed in [flLlldjl . This function, 
V(q,T), can be interpreted as the free-energy cost paid 
to keep two copies of the system at a fixed overlap q, 

V(q,T)=F q (T)-F (T); (5) 

here F q (T) is the free energy of two copies of the system 
when constrained to have fixed overlap q, and Fq(T) = 
2F(T) is the free energy of two independent copies of the 
system. 

Its qualitative behavior is shown in Fig.Q]for the 1RSB 
mean field p-spin spherical model: for T > Td it is a 
convex function of q with only one minimum at q = 0. 
At the dynamical transition temperature Td a secondary 
minimum starts to develop at finite q. On lowering the 
temperature below X^, the value of V at the minimum de- 
creases and vanishes at the thermodynamical transition 
temperature Tk- 

From the potential V(q,T) one can extract informa- 
tion about the complexity E(T) and the barrier height 
£{T). Indeed, it is well known that for T K <T <T d the 
phase space of 1RSB models is disconnected in an expo- 
nentially high number of metastable states. The Gibbs 
equilibrium state is a superposition of a subset of these 
states ( "equilibrium states" ) having a defined self-overlap 
q(T); however, the probability of finding - at equilibrium 
- two independent copies of the system in the same state 
is zero Different states have zero overlap, we there- 
fore expect that the stable phase of the two copies of the 



system - i.e., the one for which V{q) is minimum - is at 
q = 0. This is indeed the case as one can see from Fig.^ 
The secondary minimum at q ^= can then be interpreted 
as a metastable state for the two coupled systems, that 
corresponds to the situation where both systems are in 
the same state with self overlap q{T). Thus, the value 
qmin(T) where V(q, T) has a secondary minimum can be 
interpreted as the self-overlap of the equilibrium states 
at temperature T. 

The free energy of the system for Tk < T < Td can be 
written as 

F(T) = f(T) - TE(T) , (6) 

where f(T) is the free energy of a single equilibrium 
state at temperature T and S(T) is the complexity, i.e. 
the logarithm of the number of equilibrium states. The 
free energy of two independent copies of the system is 
F Q (T) = 2F(T) = 2/(T)-2TE(T), while the free energy 
of two copies constrained to be in the same equilibrium 
state is given by 

F qmin (T) = 2f(T)-Ti:(T) . (7) 

Thus, 

V{q min , T) = F qmin (T) - F (T) = TS(T) , (8) 

and the equilibrium complexity S(T) can be deduced 
from the function V(q, T). 

The difference between the value of V at the maximum 
and the value of V at the minimum can be interpreted 
as the height of a "barrier" that the two coupled systems 
have to overcome to escape from the situation where they 
are constrained to be in the same state. Thus, we can de- 
fine the "barrier height" £{T) = V{q maXl T)-V{q mm ,T). 
Note that a system-dependent proportionality factor is 
needed in order to account for the cooperativity of the 
process of escaping from a state: indeed, V(q) is the free 
energy per spin, while an unknown number of spin can 
be involved in the escaping process. Therefore, the £ (T) 
defined above is an estimate of the barrier height up to 
an unknown (system-dependent) proportionality factor. 

B. Temperatures 

The thermodynamical transition temperature Tk is 
defined as the temperature where the complexity van- 
ishes: T,(T K ) = 0. Then, at T K the value of V at the 
secondary minimum becomes equal to zero (see Fig. ^| . 
The dynamical transition temperature Td is the temper- 
ature at which the metastable minimum first appears. 
We now provide a definition of "glass transition temper- 
ature" T g . As we discussed at the beginning of this sec- 
tion, using the Adam-Gibbs relation, the usual definition 
of glass transition temperature turns out to be 

^ C (9) 
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The value of the constant C is arbitrary. Taking into ac- 
count the fact that in the considered models £ is defined 
up to a proportionality factor, we can fix the value of the 
constant in order to obtain reasonable (with respect to 
experiments) values for the different quantities we want 
to study, fragility in particular. Different choices of the 
constant change only quantitatively the results, while the 
qualitative picture stays the same. 



C. Complexity, barrier heights and fragility 

Given a definition of T g , the complexity at T g is sim- 
ply E(T fl ) and the barrier height £{T g ): clearly, these two 
quantities are related by Eq. |5| Knowing the complex- 
ity as a function of the temperature, we can define the 
fragility as: 



m(T g ) = \ + T t 



(10) 



The latter definition is very useful in a mean field context 
as - once a definition of T g has been chosen - it involves 
only the complexity, that is a well-defined quantity in 
mean field models. It is equivalent to the usual Angell 
definition of fragility if rjoo =const, and the Adam-Gibbs 
relation is assumed to be valid lj. This definition of 
fragility has been shown to be related to the one usually 
considered in experiments in Ref. 0- 



The quantity fk(T) (in the low-fc limit) can be iden- 
tified in the considered models with the self-overlap of 
the states: this identification come from the observation 
that both quantities represent the plateau of a relevant 
correlation function. Thus, we will define 



a(T g ) = 1 - q(T g ) 



(13) 



where q{T g ) is the self-overlap of the equilibrium states 
at T g , i.e., the value of q where V(q, T) has the secondary 
minimum at T — T g (see Fig. 0. 

As the self-overlap of the states is related to their vol- 
ume in phase space (high overlap corresponding to small 
states), a small value of a corresponds to small- volume 
states, while a big value of a corresponds to large-volume 
states. In this sense, a(T g ) will be called "volume of 
the equilibrium states at T g " . Note that a similar iden- 
tification has been discussed in Ref. @: indeed, from 
Eq. 7 of Ref. @ (note that due to a typing error the 
power —1 has to be disregarded) one can see that a is 
related to the curvatures of the minima of the potential 
(in the harmonic approximation), and that small curva- 
tures (large volume) correspond to large a, while high 
curvatures (small volume) correspond to small a. This 
is consistent with the equivalence of the definition of a 
given in Ref. @ and the one adopted here. 



E. Summary of the definitions 



D. Volume of the states 



As we discussed in the Introduction, in Ref. 6] fragility 
has been shown to be correlated with an index related to 
the volume of the states populated at equilibrium around 
T g . More precisely, in [y] this index has been defined as 



a{T g ) = lim 



d[fk(T)]- 



fc^O 



d(T/Tg 



(11) 



T=0 



where fk(T) is the non-ergodicity factor extracted from 
the dynamic structure factor S(k 7 to) at a given wave vec- 
tor k. From Fig. 2 of Ref. 0, we see that the possibility 
of classifying the considered systems in term of a - given 
by Eq. ^2 - relies on the observation that the curves of 
fk as a function of T/T g for different systems do not 
intersect (the same observation, that holds for log r\ (T) 
as a function of T g /T, is the basis of the definition of 
fragility). Therefore, the index a defined in [£| can be re- 
placed by other equivalent - by equivalent we mean pos- 
itively correlated - definitions (like the definition of F l j 2 
as a "fragility index" @). An useful equivalent definition 
of a is 



a{T g 



= lim 



1 - fk(T g ) 



(12) 



As one can easily check observing Fig. 2 of Ref. Q , this 
definition is equivalent to Eq. ^2 if the curves fk{T) do 
not intersect. 



To conclude this section, we give a short summary of 
all the definition we discussed. We will call q m in(T) the 
value of q where V(q, T) has the secondary minimum, and 
qmax(T) the value of q where V(q,T) has a maximum. 
Then, we define: 



E(T) 

£{T) 

T K 

T g 

T d 

miT g ) 

a{Tg) 



V(q, 

V(q 

S(T K 



(T),T)/T 

tn ax it 

= 



>(T),T) 



£(T S ) 



C 



Qmax (T d ) = 



(Td) 



T 9 E(T g ) 
qmin \T g ) 



The constant C will be chosen in order for the fragility 
to be in the experimentally observed range. 



III. SPHERICAL P-SPIN MODEL 

In this section, we will compute explicitly all the pre- 
viously defined quantities in the spherical p-spin model. 
The model is defined by the Hamiltonian 



Ju 



(ii,— i»j>) 



(14) 



where <7j are real variables subject to a spherical con- 
straint J^CT? = N, and J^...^ are quenched ran- 
dom Gaussian variables with zero mean and variance 
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FIG. 2: Thermodynamic transition temperature Tk, glass 
transition temperature T a and dynamical transition temper- 
ature Td for the p-spin spherical model as a function of p — 2. 



p\/(2N p ~ 1 ). This simple model has been successfully 
used for studies of the glass transition 0, 0|- It is a 
"Gaussian-like" model, in the sense that its complexity - 
even if the distribution of states in not exactly Gaussian 
- is known to be a concave function of the temperature, 
that vanish at Tk and assumes its maximum at Td, with- 
out any inflection point in between |13| . 

The expression for V(q,T) in the p-spin spherical 
model has been computed in Refs. |3, Il4j . However, 
a simplified expression can be used when the value of 
V(q,T) on its stationary points is considered (see Ap- 
pendix [SJ : 



V(q,T) 



T 



log(l - q) 



Tq 
2 



(15) 



This function can be shown to coincide with the cor- 
rect V(q,T) on each stationary point of V(q,T). As we 
are interested only in the value of V(q, T) on its station- 
ary points, the use of the correct V(q, T) calculated in 
Refs. 0, uM or 01 the one given by Eq. ^] will give ex- 
actly the same result. 

Note that, while the model is defined only for integer p, 
Eq. E| holds also for real p; we will therefore discuss the 
behavior of the different quantities for any real p > 2. In 
particular, the p — > 2 limit is interesting being related to 
a diverging fragility (Td — ► Tk) and to the discontinuous 
1RSB transition becoming a continuous one. 



A. Temperatures 

From Eq. 1151 we can compute the three temperatures 
Tk, T g and Td as functions of p. Their behavior is re- 
ported in Fig. |3 We immediately note that, for p ~ 2, 
the difference between Tk and T g is very small, therefore 
the system is very fragile; moreover, for p — > oo the Kauz- 
mann temperature approaches zero (as 1 / ' ^/\ogp) , while 



FIG. 3: The complexity E(T 9 )/£(T) as a function of T a /T for 
the p-spin spherical model at different values of p. Fragility 
is the slope of the curves in T a /T = 1. The system becomes 
stronger on increasing p. 



the glass transition temperature remains finite. The sys- 
tem therefore becomes stronger and stronger on increas- 
ing P- 



B. Complexity and fragility 

The same observation can be made more quantitative 
by considering an " Angell plot" for the complexity Q : in 
Fig. we show the complexity S(T) as a function of the 
temperature, for different values of p. The choice of the 
particular scaling that appears in Fig. |3| has been made 
in order to make a close correspondence with Fig. 2 of 
Ref. 0- We see that the curves for different values of 
p are ordered from bottom to top. The same behavior 
is observed in experimental systems of different fragility. 
Indeed, the index of fragility defined in Eq. ^|is exactly 
one plus the slope of the curves in T g /T = 1 (see Fig. EJ : 



m(T g ) = 1 + T 



= 1 



d[s(r fl )/s(r)] 



d[T g /T] 



(16) 



T=T„ 



The fragility index m is shown in Fig. 0] as a function 
of p. We see that it is a decreasing function of p. Its 
values are in the range observed for experimental system 
due to our (arbitrary) choice of the constant C appearing 
in Eq. El C = 0.1. In Fig. |E(T S ) is also reported as a 
function of p. We see that it is an increasing function 
of p, that diverge as logp for p — > oo: thus, the number 
of states in this system is a decreasing function of the 
fragility, at variance with what is currently believed (for 
a review, see Ref. 0). We will discuss this point in detail 
in section Ivl 
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FIG. 4: Fragility m(T 9 ), configurational entropy E(T S ), "vol- 
ume" of the equilibrium states a(T g ) and barrier height S(T g ) 
for the p-spin spherical model as a function of p — 2. 



FIG. 5: Thermodynamic transition temperature Tk, glass 
transition temperature T g and dynamical transition temper- 
ature T d for the p-spin Ising model as a function of p — 2. 



C. Barrier heights and volume of the states 

In Fig. ^the barrier height £(T g ) is also reported as a 
function of p, together with the index a(T g ) = 1 — q(T g ) 
that we called "volume" of the equilibrium states at T g . 
We observe that in this model the states become smaller 
on increasing p, while the barriers separating them in- 
crease. In section [V] we will discuss this behavior try- 
ing to deduce a geometric description of the evolution of 
the phase space of this model at different p, and relate 
fragility to geometric properties of the phase space. 



IV. ISING P-SPIN MODEL 

The Ising p-spin model is anoth er p opular model for 
the study of the glass transition [15|,[l6|]. Its Hamiltonian 
is given by Eg. 1141 where the variables are Ising spins, 
cr, = ±1, and the spherical constraint is absent. For the 
Ising p-spin model, the two-replica potential V{q,T) is 
given by 

V(q,T) = P^-q P + P-q?- 1 

J Vz cosh(Az) logcosh(Az) v ' 
J Vz cosh(Az) ' 

where Vz = exp(-z 2 /2) dz, and A 2 = /3 2 §q p_1 . 

The Ising p-spin model is also a "Gaussian-like" model, 
like the spherical one. However, the total number of 
states in the Ising p-spin model cannot be greater than 
2 N (the total number of configurations), and hence 
E(T) < log 2, while in the spherical model S(T ff ) diverge 
as logp for p — > oo, as previously discussed. 



A. Temperatures 

The first consequence of this difference is observed 
when studying the transition temperatures as a function 
of p (see Fig. |SJ). Indeed, as in the spherical model, we 
have Tk ~ T g for p ~ 2, and T g ^> Tk for p — ► oo. 
But, in this model, Tk tends to a finite value at large 
p, while T g and Td diverge. This behavior can be under- 
stood recalling that for a "Gaussian-like" model we have 
Tk ~ l/VS*) beingthe total number of states, i.e., 
the maximum of E(T) |l|. 

B. Complexity and geometric properties of the 
phase space 

The "Angell plot" for the complexity of the Ising p- 
spin model looks very similar to the one of the spherical 
model (see Fig. |3J) and is not reported here. 

Having fixed an appropriate value for the constant C 
in Eq. |U(C = 0.02, different from the value chosen in the 
previous case) , the behavior of the fragility as a function 
of p is also very similar to the one of the spherical model. 
The same behavior is found for the other quantities under 
study, as one can deduce from a comparison of Fig. Eland 
Fig. 0] the main difference being the discussed behavior 
of S(T 9 ) at large p. 

C. Vibrational properties and volume of the states 

Another relevant difference between the spherical and 
the Ising model is that, in the latter, harmonic vibrations 
are not present (the variables being discrete): we have 
q(T) — > 1 exponentially for T ~ 0, and the definition of 
a via Eq. gives a = for all p. However, the defini- 
tion given in Eq. 1131 and used in our calculations gives a 
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FIG. 6: Fragility m(T 9 ), configurational entropy E(T S ), "vol- 
ume" of the equilibrium states a(T g ) and barrier height S(T g ) 
for the p-spin Ising model as a function of p — 2. 
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FIG. 7: Fragility versus a for the two investigated models. 
The curve is very similar for the two models, and is consistent 
with the linear correlation found in at least for a < 0.4. 



reasonable result also in absence of harmonic vibrations. 



V. CORRELATIONS BETWEEN DIFFERENT 
PROPERTIES OF THE PHASE SPACE 

In this section we will examine the correlations in the 
quantities under study, trying to relate fragility to the 
phase space geometry. We will compare our results with 
the general consideration that we made in Ref. Q], and 
with the experimental results of Ref. |(j . 



A. Fragility and volume of the states 

In Ref. |y| it has been established that fragility is pos- 
itively correlated with the index a defined in section ITT1 
In other words, fragile systems have large basins while 
strong systems have small basins. In Fig. [7| we plot the 
fragility m as a function of a parametrically in p for the 
investigated systems. The curve m(a) is very similar for 
the two models - remember that the only adjustable pa- 
rameter is the constant C in Eq. [5] By comparison with 
Fig. 3 of Ref. , we conclude that the model has a be- 
havior similar to the one of real systems. Surprisingly, 
also the linear correlation between m and a is reproduced 
for a < 0.4. Thus, mean field p-spin models are able to 
describe the relation between fragility and the volume of 
the basins visited around T g found in Ref. [||. 



However, in the model considered here the behavior is 
exactly the opposite. In Fig. |H1 we report S(T S ) as a 
function of the fragility: we see that the total number of 
states is a decreasing function of the fragility, at variance 
of what is currently believed. To discuss this point, we 
have to refer to Q: there, we discussed the possibility 
of correlating fragility with the total number of states 
for general models of S(T), and assuming the validity of 
the Adam-Gibbs relation, Eq. ^ We concluded that the 
knowledge of the distribution of states is not enough to 
determine the fragility. Indeed, the relevant parameter 
was identified, for a general "Gaussian-like" distribution 
of states, as 



D = 



£(T g ) 
T K Z(T g ) 



(18) 



Note that in Eq. ^]we have to calculate £ at T = T g be- 
cause in the considered models the barrier height £ is a 
T-dependent quantity, while in the Adam-Gibbs relation 
it is usually assumed to be a constant (see Eq. Q). How- 
ever, the Adam-Gibbs relation has been tested around 
T g , therefore, to a good approximation, we can fix £ to 
be a constant equal to its T = T g value. The parameter 
D is inversely proportional to the fragility m: therefore 
m ~ T,/£. Thus, fragility is not simply correlated to 
the total number of states: if the "barrier heights" grow 
faster than the total number of states, fragility can be a 
decreasing function of S. We will now show that this is 
indeed the case in the considered models. 



B. Fragility and total number of states 

It is usually believed that fragile systems have a larger 
number of states than strong ones, even if the total num- 
ber of states is not an experimentally accessible quantit; 
and numerical simulations give contradictory results 



rtity 

m 



C. Barrier heights, total number of states and 
fragility 

From Fig. 0] and Fig. we see that the barrier height 
is indeed an increasing function of p in the considered 
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FIG. 8: Total number of states (represented by the complexity 
at Tg) as a function of the fragility m: an inverse correlation is 
found between these quantities, at variance of what is naively 
expected from the Adam-Gibbs relation. 
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FIG. 9: Sketch of the evolution of the p-spin free energy by 
varying p: at small p there is a small number of states of large 
volume separated by low barriers; at high p there is a large 
number of states of small volume separated by high barriers. 
The height of the barriers increase faster than the number of 
states: thus, fragility is a decreasing function of p. 



models. Using Eq. Eq. EDcan be written as 

D=C^. (19) 

Therefore, from Fig.|21and Fig.[S]we see that D is indeed 
an increasing function of p that diverge for p — > oo, as 
the ratio T 9 /Tk increase on increasing p for both models. 
Thus, we can conclude that in the considered models the 
height of the barriers (in units of Tk) increases faster 
than the total number of states. This explains why one 
observe an inverse correlation between fragility and the 
total number of states, as discussed above and in Ref. Q. 

D. A geometric picture of the phase space 

Collecting all the information that we obtained in the 
previous sections, we can propose a geometric picture of 
the variation with p of the p-spin model free energy land- 
scape. Indeed, on increasing p: 

i) The total number of states increases. 

ii) The volume of the states decreases (a decreases). 
Hi) The height of the barriers between states increases. 
Thus, we get the picture of a landscape where, on in- 
creasing p, a great number of small states with very high 
curvatures and separated by very high barriers appear: a 
sketch of this evolution is given in Fig. [5| The behavior 
of the fragility in this situation is related to the behavior 
of £/£, the ratio between number of states and height 
of the barriers between them: in these models, it turns 
out that £ increase faster than E, and the fragility is a 
decreasing function of p. 

This behavior is consistent with the fact that fragility 
turns out to be positively correlated with the "volume" of 
the states as measured by a. Indeed, if, on the contrary, 
the barrier height grew slower than the total number of 



states (equivalently, if m would be positively correlated 
with the total number of states), there should be also 
an inverse correlation between m and a, in disagreement 
with what is experimentally observed. 

In the p — > 2 limit, where the fragility becomes infinite, 
the second derivative with respect to q of the potential 
V(q,T) calculated in q = and T = Tk — Td vanishes 
(see Fig. ^| and the so-called spin glass susceptibility di- 
verges at the critical temperature. In other words when 
the fragility becomes infinite soft modes appear at the 
critical temperature supporting the previously presented 
physical picture. 

Note that the outlined picture is valid for "Gaussian- 
like" models, i.e., models where the complexity is a con- 
cave function of the temperature that vanish at Tk with- 
out any inflection point. These models seem to describe 
correctly the distribution of basins in real systems only 
for relatively high fragilities. The behavior of the com- 
plexity (or configurational entropy, or excess entropy) as 
a function of temperature for very strong systems is still 
an open problem; our discussion may not apply to these 
systems. 

VI. CONCLUSIONS 

From the investigation of two mean field p-spin models, 
we proposed a picture for the evolution of the free-energy 
landscape from fragile liquids to strong ones. This pic- 
ture accounts for the recently observed correlation be- 
tween fragility of a liquid and vibrational properties of 
the corresponding glass. The main prediction of our anal- 
ysis is that the total number of states and the Adam- 
Gibbs parameter £ should both be decreasing functions 
of the fragility. Unfortunately, existing data are not suffi- 
cient to strictly test this prediction; the excess entropy is 
available only for few experimental systems, and numer- 
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ical simulations are performed in a temperature range 
where the fragility of the investigated systems is approx- 
imately the same. We hope that this predictions can be 
tested in the future. 

We thank Stefano Mossa for a careful reading of the 
manuscript, and Cristiano De Michele, Tullio Scopigno, 
Francesco Sciortino and Luca Angelani for useful discus- 
sions. 



APPENDIX A: THE TWO-REPLICA 
POTENTIAL 

The two-replica potential is defined in Ref. as 



V(q) = —F(T) - - J ^-^ffT log Z(a, q) 
Z(a,q) = [ dre-^5( q - q (a,T)) , 



(Al) 



where F(T) is the equilibrium free energy, and q(a,r) is 
the overlap function. The following expression is then 
derived (la ]: 

V(q) = —F(T) 

_ Um i im 2-JL ( f dae-WMzfaq)*"- 1 )' 



(A2) 



The last integral can be rewritten as 



J dae-P H ^Z{a,q) m - l X = 



„ n m 

/ da aa e-^™ H ^H Y[S(q-q(<T al ,a aa )) , 



(A3) 



where a — 1, • • • , n, a — 1, • • • , m. This is exactly the 
expression of the nm times replicated equilibrium par- 
tition function, with the additional constraint given by 
the 5-functions. Using standard manipulations |l2(, it is 
rewritten as 

/n m 
dQ aa ,bi3 e w/(Q) J] J] S(q- Q ahaa ) 

a=la= \ (A4) 



Z 



Thus, evaluating the integral at the saddle point, we get 



V(q) = -F(T) - lim lim ~^f(Q) . (A5) 
n^o m— >i n am 

The matrix Q is defined by the following conditions: 

i) the elements on the diagonal are equal to 1; 

ii) the elements Q a iaa, a > 1, are equal to q; 

Hi) all the other elements are determined by the maxi- 
mization of f(Q). 



As usual, one needs a parametrization of the matrix Q in 
order to perform the analytic continuation to non-integer 
n and m. A possible ansatz is 9] (in the example, n = 3, 
m = 4): 

\ 



Q 





V 



Within this ansatz, and using the relation 





(l-r) m - 2 [l-2r + rm-(m-l)q 2 } , 



T 
2 



log(l - r) + 



r — q 
1-r 



(A6) 



, (A7) 



where r(q) is determined by d r V = 0. Now, it is easy 
to check that the condition dV/dq = d q V = 0, together 
with d r V = 0, is satisfied if q = r. Thus, when V(q) is 
stationary, r(q) = q and the potential V(q) reduces to 
the one given by Eq. 

The fact that when dV/dq = the matrix Q reduces to 
the usual 1RSB overlap matrix (let us call it Q) is general: 
indeed, the condition dV/dq from Eci. lA5l is equivalent to 



df(Q) 
dQ 







(A8) 



This means that the function f(Q) must be stationary 
with respect to all the elements of Q if dV/dq — 0, and 
we know that the 1RSB matrix Q provides a solution to 
this condition. As a final remark, we note that if Q = Q, 
we have 



f(Q) 



iRSB{m,q) 



(A9) 



where 4>irsb is the usual 1RSB free energy. Substituting 
this expression in Eq. IA5I one obtains 



V(q) - -F(T) + 



lim d n 

re->l 



mtpi R sB{m 



lim d m <j) 1R sB(m,q) 

m— *1 



(A10) 



using the relation 4>iRSB( r n = 1) = F(T) that holds 
above Tk- Therefore, on its stationary points, V(q) is 
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given (at the 1RSB level) by this simple expression, that for the computation of V(q). However, we don't account 
can be easily calculated in several models. Note that, as for them in this paper, 
discussed in Ref. [T3 [. full RSB effects can be important 
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